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The measurements of flow harmonics V2-vq using the event plane and two particle correlations methods in broad 
Pt, V and centrality ranges using the ATLAS detector at LHC are presented. ATLAS recorded about 9 /ib~^ 
of lead-lead collision data in the 2010 heavy ion run. The full azimuthal acceptance of the ATLAS detector 
in ±2.5 units of pseudorapidity for charged hadrons and the large amount of data allows for a detailed study 
of the flow harmonics. The pj-, centrality and r] ranges where the two methods give consistent Vn and where 
they disagree are discussed. It is shown that the ridge as well as the so called "mach-cone" seen in two particle 
correlations are largely accounted for by the collective flow. Some scaling relations in the Pt dependence of the 
Vn are also discussed. 



1. Introduction 

A central goal of the heavy ion program at the LHC is to understand the properties of the hot and dense 
matter produced in relativistic heavy ion collisions, commonly termed as the quark gluon plasma. One of the 
ways to study its properties is to perform high accuracy measurements of the azimuthal yields of the produced 
particles. Due to anisotropic pressure gradients, the "fireball" produced in the collisions expands differently in 
different directions with more particles emitted along the larger gradients. This anisotropy in the particle yields 
can be expressed as a Fourier series: 



dN 



(xl + 2Y,VnCOsini(j)-^n)) (1) 



Where (j) is the azimuthal angle of the produced particles, Vn is the magnitude of the n*'* order harmonic 
flow and is the orientation of the corresponding harmonic plane The superscript 'a' is used to indicate 
particle species, pt bin, centrality, etc. 

The elliptic flow V2 has been studied extensively as it usually dominates other coefficients, being strongly 
influenced by the elliptic geometry of the overlap region between the colliding nuclei. However, the initial 
geometry of the produced fireball may have angular moments of several orders: elliptic, triangular, etc due to 
fluctuations in the spatial positions of the nucleons in the colliding nuclei. These higher order moments in the 
initial geometry can give rise to measurable higher order flow (n > 3). 
" ■ A better understanding of the higher order Vn can explain the origin of the ridge, an elongated structure 

along A77 at Acj) ~ f^l and the so-called "mach-cone" , a double hump structure on the away side 3 seen in 
two particle correlations. These were initially interpreted as response of the medium to the energy deposited by 
the quenched jets. However, recent studies [J] have shown that higher order flow harmonics make a significant 
contribution to these structures. The Vn can also be used to obtain information about the initial geometry and 
transport properties of the medium, for instance the viscosity to entropy ratio: ry/s 

The Vn can be measured by the event plane (EP) method where one determines the orientation of the harmonic 
planes \l/„ and then, by measuring the yield of the particles about these planes, obtains the u„. They can also 
be measured by the two particle correlation (2PC) method where one measures the relative yield of associated 
(or partner) particles with respect to trigger particles in A(j> = 4>trigger — 4'partner and Arj ~ rjtrigger — Vpartner ■ 
The 2PC correlation function can be expanded in a Fourier series in Acj) as: 



C{A^) oil + 2Y,Vn,n{PT,PT)cos{nA^) (2) 

n=l 

Where the superscripts a and b label the trigger and partner particles. If the only contribution to correlations 
comes from collective flow, then the Fourier coefficients Vn,n are equal to the product of the individual single 
particle Vn [1]: 



(3) 
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Using EqISl one can obtain the w„ from the 2PC. The above relation is violated if non-flow effects, like jets, 
are present. Thus before extracting the flow harmonics from the correlations, one must ensure that most (if not 
all) of the non-flow effects are removed. 

All the results presented in this proceeding are for inclusive charged hadrons reconstructed using the ATLAS 
inner detector covering the pseudorapidity range \rj\ < 2.5. For the event plane analysis, the harmonic 
planes are determined using the Q-vector method The Q-vectors are measured using the ATLAS forward 
calorimeter (FCal) which lies at G (3.3,4.8). Details of the results presented here can be found in [lol [Tl|. 



2. Event plane analysis 

The measured vary around the true due to statistical fluctuations (H. In order to correct for this, 

the Vn are obtained by first calculating = (cosn(0 — V&n)) and then dividing the v^'^ by the reaction plane 
resolution 



Vn = v°^y Resolution{'ii n) (4) 

The FCal reaction plane resolutions for n=2-6 are shown in FiglTJ For different harmonics, the resolutions 
are cut off at different centralities depending on where they become small as compared to the systematic and 
statistical uncertainties. 



■J _ /I rZ./lS Preliminary |Ldt = 8nb' ^n=2 _ 




Central ity [%] 



Figure 1: FCal reaction plane resolutions for harmonics 2 to 6. Bands indicate systematic uncertainties. 



In FigHlthe w„ are plotted as a function of centrality from EP method in 5% centrality bins plus a 0-1% super- 
central bin. Each panel is a different px range. It is seen that V2 has a much stronger centrality dependence 
as compared to the other harmonics. For non-central events, V2 becomes much larger than the other v„ due to 
the elliptical geometry of the fireball. The other Vn, mainly driven by event-by-event fluctuations in the initial 
geometry have a weaker centrality dependence. This centrality dependence of the collision geometry is also the 
reason why the reaction plane resolution for ^['2 shown in FigIT] increases from central to mid central events, 
while the resolution for the other do not show such a strong centrality dependence. Also note that in most 
central events (top 5% or top 1%) and at sufficiently high px, V3, V4 and even ^5 can become larger than or 
comparable to V2- 

Figure [3] shows Vn as a function of pt- For all centralities the Vn have a similar trend, they increase to reach a 
maximum between 3-4 GeV and then decrease. In the 0-5% central events, V3 and U4 become larger than V2 at 
intermediate px- A scaling, previously observed at RHIC 12| is seen in the px dependence of ?;„: UtV" = knV^^^ 
where fc„ are only weakly dependent on px- In Fig|3]the values of fc„ are plotted as a function of px for different 
centralities. The fc„, though weakly dependent on px, have a strong centrality dependence. 

In Figl5] the i] dependence of the v„ are shown for 2 < < 3 GeV. In this px interval, V2 decreases by at 
most 5% from 77 = Oto?7 = 2.5, for n > 2 the decrease is slightly larger. In general, there is a weak dependence 
of the differential Vn on rj in central and mid-central collisions. For V2, this weak 77 dependence is demonstrated 
over a large px and centrality range in Fig|51 
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Figure 2: u„ as a function of centrality is 5% centrality bins plus a 0-1% most central bin (the right most point) for four 
Pt bins. The shaded bands indicate systematic uncertainties. 
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Figure 3: Vn vs. pr for several centrality classes. The shaded bands indicate the systematic uncertainties. Note that the 
ve is only measured for 0-40%. 
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Figure 5: dependence of u„ for pr G (2, 3)GeV. Each panel is a different centrality bin. 



Eq|3] stated that for pure flow, the 2PC is the product of the single particle Vn- This relation is only true 
if there is no r/ dependence in the u„ (or approximately true if the r] dependence is weak). Otherwise Eq|3] would 
have to be replaced with a convoluted integral of the rj dependent Vn weighted by the 77 dependent spectra. The 
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Figure 6: 77 dependence of V2 for 0.5 < pr < 20 GeV in five pr intervals and 10% centrality intervals. Error bars show 
statistical and systematic uncertainties added in quadrature. 



3. Two particle correlations 

Figure [7] shows the 2PC in Arj-Acp space measured in the momentum range 2 < p^,p^ < 3 GeV for several 
centrality classes. The near side jet peak reported in 0] is clearly visible at (At;, At/)) ~ (0, 0) across all 
centralities and is more prominent in peripheral events. The near side ridge has a weak Arj dependence in 
central and mid-central collisions. The strength of the ridge increases from central to mid-central events and 
then decreases for peripheral events, completely disappearing in the 80-90% centrality class, leaving behind only 
the near-side jet peak. In the most central events (first and second panels) one can see a double hump structure 
on the away side which transforms into a ridge-like structure in the mid-central events. In peripheral events, it 
is replaced by the away side jet. 

The presence of the near side jet biases the values of the flow harmonics obtained from the two particle 
correlations. This bias can be minimized by requiring a |A?7| > 2.0 cut between the trigger and the associated 
particle in the correlations. While this cut removes the near side jet, part of the away side jet is still left. This 
is because the away side jet is not localized in Arj and cannot be removed by a cut on IA77I. Thus the measured 
Vn,n still retain non-flow bias. 

Figure |8] explains the procedure used for obtaining the w„ from the 2PC. The panel a) shows the two- 
dimensional correlation function in Arj-Acj) for 2 < < 3 GeV. Such correlations, where both trigger and 
partner have same pT range, are termed as fixed-pT correlations. The panel b) shows the the one-dimensional 
A(j) correlation function (black points) in Ac/) for 2 < |Ar/| < 5, overlaid with contributions from individual 
Fourier components and their sum. As mentioned before, the IA77I > 2.0 cut is imposed to remove the near 
side jet. The Fourier coefficients Vn,n obtained from the one-dimensional correlation are shown in panel c) as a 
function of n. The bars and bands indicate statistical and systematic uncertainties respectively. As the trigger 
and partner px ranges are chosen to be the same, Eql3] becomes: 
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Figure 7: Two particle A77 — correlations for 2 < PtiPt < 3GeV. Each panel is a different centrality bin. The 
peripheral events have the near side jet peak truncated. 



VuApt^Pt) = «'(Pt) (5) 

Using the above relation one calculates the values of w„ and its uncertainties from The w„ values are 

plotted as a function of n in the panel d) of Fig 15] In panels c) and d) of Fig 15] the u„^„ and Vn values are plotted 
untill n = 15. However, the analysis is limited to n < 6 because for higher n the systematic and statistical 
uncertainties are large compared to the extracted value of Vn,n- 

As mentioned before, a |A?7| > 2.0 cut was used to remove the effects of the near side jet. The choice of this 
particular value of the cut can be justified in the following way. Figure [5] shows the and w„ obtained for 
2-3 GeV fixed-pT correlation as a function of |Ary|. The Vnn and Vn show a clear peak at jAyyl ~ due to the 
near side jet. For larger values of jAryj, the Vn,n values drop as the influence of the near side jet weakens. For 
jAryl > 2, the fractional change is very small. 

The two-particle correlation analysis used the fixed-pT correlations and Eql5] to obtain the w„ values from 
the Vn.n- The Vn values obtained from fixed-p^ can be cross-checked using correlations where the trigger and 
partner have different and values. Such correlations are termed as mixed-px correlations. The Vnip^x) '^^ 
the partner can be obtained as: 

VniPx) = Vn,n{PT:PT)/'"n{PT) (6) 

Equation [6] can be checked by measuring the same for different w„(p^). This is illustrated in FiglTO] In 

the top panels the V2 — f 5 for four different pr obtained from fixed-pr correlations bins are plotted as a function 
of I A77I for the 10-20% centrality events. For these four values of the trigger pT, mixed-pT correlations are done 
with the associated particle p^ S (1.4.1.6) GeV. The bottom panels of FigflUl show the Vn{Px) obtained from 
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Figure 8: The steps involved in the extraction of the ii„ for 2-3 GeV fixed-pT correlation. Panel a): the two-dimentional 
correlation function (shown for \/S.rj\ < 4.75 to reduce fluctuations near the edge). Panel b): The one-dimentional A(j!) 
correlation function for 2 < \Arj\ < 5. Panel c): v„^n vs. n. Panel d): Vn vs. n. The lower sub panels in c) and d) also 
show the v„^n and ii„ respectively for n > 7 but in a linear scale. 
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Figure 9: Vn,n and v„ as a function of A?7 for 2-3 GeV fixed-pT correlation for events in 0-5% centrality class. 



the mixed-pT correlations using EqlHl At low IA77I different triggers give different values for Vnip^) sfiowing that 
the scaling relation Eql3](or Eql6|) breaks down due to the effects of the near side jet. For \ Ar]\ > 1.0 however, 
all trigger-partner combinations give nearly identical values of Vnip^) even though the trigger Vn values vary 
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over a large range. This validates EqI3] and EqlBl and shows that indeed the Fourier components in the 2PCs 
at large jAryl are due to collective flow. Note that the mixed-pT correlations are only used as a cross-check to 
qualitatively check that Eql3] works, they are not used to determine any systematic errors for the 2PC Vn- 
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Figure 10: Top panels: v„{pt) ~ \/\'Vn,n{PTtPT)\ {n = 2 — 5 from left to right) vs. |A7;| for four fixed-pT correlations 
{Pt G 0.5-1, 1-2, 2-3 and 3-4 GeV ranges). Bottom panels: VniPx) = 'fn,n{pT)PT)/^ri (p^) vs. lAr;! forp^ G (1.4, 1.6)GeV. 



4. Comparison of Vn obtained from 2PC and EP methods 

Figure [TT] shows the centrality dependence of the u„ obtained from the two methods. For clarity only the 
statistical errors are shown. From the ratio plots in the bottom panels one can see that the values for V2, V3 
and W4 agree within 5% for central and mid central events. For vc, and vq they agree within 10% and 15% 
respectively and are consistent within the systematic uncertainties quoted in 10]. 

In Fig |12l the px dependence of V2 obtained from two-particle fixed-pT correlation method in several A77 slices 
are compared to the EP values. 

The 2PC values measured for jAryj S (0,0.5) for all pr are signicantly higher than the EP results due to 
the bias introduced by the near side jet. However, the curves measured with |A?7| > 2.0 agree much better 
below pt=4.0 GeV. For higher values the away side jet presumably affects the 2PC results which explains their 
deviation from the EP Vn even for jAr^l > 2.0. 

For central and mid central events and < 4 GeV, the two particle correlation with IA77I > 2 gap and EP 
method give consistent values of Vn ■ This indicates that the structures in two particle correlations at large Arj 
including the ridge and cone can be entirely accounted for by collective flow. To demonstrate this clearly, the 
EP Vn values are used to reconstruct the 2PC which are then compared to the measured 2PC. The reconstructed 
Vn^n are calculated according to Eql3]as the product of the Vn obtained using the EP method: 



<r = X V^ (7) 

Since Vi measurements were not made using the EP method, the wi^i component of the reconstructed corre- 
lation function is taken to be the same as that of the measured correlation function. The overall normalization 
of the reconstructed correlation function is taken to be the same as that of the measured correlation function 
{Nq^'^). The final expression for the reconstructed correlation function can be written as: 



C'^^^°(A(/.) = TV^PC / ^ ^ 2vl'l^ cos(A(/.) + 2 ^ vf/v^'' cos(nA0) J (8) 

Where, as explained before the normalization and the vi^i component are taken to be the same as that of the 
measured correlation. 
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Figure 11: Comparison of ii„ obtained from the two-particle fixed-pr correlation method (solid symbols) to the EP 
results (open symbols). The ratios between the two methods are shown in the bottom panels. Results for two pr bins 
are shown: 1 — 2GeV (left panels) and 2 — 3GeV (right panels). The dashed lines indicate a ±10% band to guide the 
eye. 
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Figure 12: Top panels: Comparison of V2 determined using the EP method and the two-particle fixed-pr correlation 
method in several A?7 slices for three centrality bins. Bottom panels: Ratios of V2 between the two-particle and EP 
results. 



In Fig ll3l four reconstructed correlation functions are compared to the corresponding measured ones. The 
2PC are well reproduced with both the near side ridge as well as the away side double hump. 



5. Summary 

To summarize, the flow harmonics V2-vq are measured by event plane method as well as two particle correlation 
method. These measurements were done to higher orders, n=6 and over a much larger pT and rj ranges than 
previously done at RHIC. An approximate scaling relationship between the w„: vU^ cx which was observed 
at RHIC was shown to hold also at LHC energies. 

For |A?7| ~ the Vn values obtained from 2PC are significantly different than the EP values, due to the effects 
of the near side jet. However, as the r] gap between the trigger and partner is increased, the 2PC f„ values 
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Figure 13: Measured correlation function data (circles) compared with that reconstructed from from two-particle 
correlations and V2 — from EP method for 0-1% centrality for several fixed-pr correlations. 



become consistent with the EP ones for central and mid-central events. 

It was demonstrated that the ridge and cone seen in two particle correlations at low and intermediate pt 
{pt < 4.0 GeV) can be entirely accomited for by the collective flow of the medium. There is no need for any 
additional physics models to account for these structures. 
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